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ABSTRACT 


The  transonic  indicial  theory  is  used  to  predict  the  effect 
of  a  gust  on  an  airfoil  at  transonic  speeds.  The  effect  of  oper¬ 
ating  two  control  surfaces  is  also  modeled  by  the  indicial 
method.  The  transonic  indicial  method  is  linear  in  a  strained 
coordinate  system  and  superposition  can  be  used.  This  allows  the 
effects  of  an  arbitrary  gust  and  control  surface  deflection  to  be 
modeled  simply  if  the  indicial  responses  for  the  gust  and  each 
control  surface  are  known.  The  computation  time  is  small  and, 
therefore,  an  optimization  technique  can  be  used  to  determine  the 
best  control  surface  deflections  to  alleviate  the  gust  loading. 


1.  Introduction 


One  of  the  more  pressing  problems  in  designing  a  modern 
aircraft  that  flies  at  transonic  speeds  is  the  prediction  of  the 
gust  response.  A  knowledge  of  the  loading  induced  by  a  gust  is 
necessary  for  the  aerodynamicist  to  predict  its  effect  on  wing 
performance  and  for  the  structural  engineer  to  design  a  suffi¬ 
ciently  strong  structure  and,  perhaps,  to  implement  active  con¬ 
trols  to  reduce  the  load.  At  present,  the  capability  to  accu¬ 
rately  predict  gust  loadinq  at  transonic  speeds  does  not  exist. 
It  is  this  problem  that  is  addressed  in  the  present  work. 

Aeroelastic  calculations  are  frequently  performed  using 
linear  theory  for  both  the  structure  and  the  aerodynamic  loads. 
This  leads  to  a  relatively  simple  means  of  estimating  flutter 
speeds,  since  linearity  permits  the  use  of  the  principle  of  su¬ 
perposition.  Linear  aerodynamic  theory  is  really  only  applicable 
to  the  subsonic  or  supersonic  speed  regimes;  the  transonic  regime 
is  essentially  nonlinear.  This  nonlinearity  is  associated  with 
shock-wave  formation,  often  accompanied  in  unsteady  flow  with 
complex  shock  oscillations.  It  is  essential,  therefore,  to  in¬ 
clude  the  nonlinearities  in  some  way  in  order  to  predict  tran¬ 
sonic  flow  phenomena.  For  the  most  general  case,  this  requires 
the  development  of  a  three-dimensional  transonic  flow  computer 
code  for  predicting  unsteady  effects  including  viscous  interac¬ 
tions  between  any  shock  waves  that  form  and  the  boundary  layer  on 
the  surface  of  the  structure.  Although  progress  is  being  made 
toward  the  development  of  such  a  code,*  this  is  only  half  the 
problem,  since  the  aerodynamic  loads  must  interact  with  a  dynamic 
structural  model  in  order  to  predict  aeroelastic  phenomena. 
Accomplishing  this  interaction  with  a  large  and  complex  unsteady 
flow  transonic  code  is  not  a  trivial  matter,  and  major  advances 
in  computer  size  and  speed  will  be  necessary  before  calculations 
of  this  sort  can  be  routinely  used  in  either  the  analysis  or  the 
design  of  aircraft. 


A  further  problem  is  that  conventional  gust  analysis  may 
impose  an  upwash  of  finite  length  into  the  free  stream.  If  this 
upwash  is  completely  arbitrary,  the  flow  is  unlikely  to  be  irro- 
tational  throughout  the  gust.  Since  all  of  the  available  tran¬ 
sonic  solution  methods  are  based  on  potential  ( irrotational ) 
theory,  '  this  presents  difficulties  in  easily  developing  a  gust 
load  prediction  method. 

If  the  gust  loading  can  be  predicted,  it  is  possible  that 
some  means  to  reduce  the  loading  can  be  developed,  perhaps  by  the 
deflection  of  a  control  surface.  Also,  due  to  the  coupling  of 
the  mean  steady  and  unsteady  loading  in  a  nonlinear  transonic 
theory,  certain  wings  may  be  more  sensitive  than  others  to  gusts. 

The  work  reported  is  concerned  with  a  preliminary  study  of 
transonic  gust  loading  and  its  alleviation.  To  retain  simplicity 
at  this  stage,  only  two-dimensional  problems  are  considered. 

In  recent  years  the  numerical  simulation  of  steady  transonic 
flow  has  reached  a  fairly  advanced  stage.  For  example,  the  two- 
dimensional  codes  based  on  the  original  algorithm  of  Murman  and 
Cole'*  are  freguently  used  in  the  aircraft  industry.  Three- 
dimensional  variations  of  these  codes  for  finite  wings  and  wing- 
body  combinations  are  also  available,  such  as  the  codes  of  Bailey 
and  Ballhaus^  (transonic  small  disturbance)  and  Jameson  and 
Caughey®  (full  potential).  While  these  codes  have  certain  defi¬ 
ciencies,  they  have  been  extremely  useful  in  aircraft  design. 

The  state  of  the  art  regarding  codes  for  unsteady  flow  is 
not  as  advanced  as  that  for  steady  transonic  flow.  The  direct 
integration  method  of  Ballhaus  and  Goorjian^  for  low  freguency 
two-dimensional  flows  is  a  reliable  tool.  There  have  been  sev¬ 
eral  extensions  of  this  theory  to  high  freguency  flows,  notably 
by  Rizzetta  and  Chin®  and  by  the  XTRANL  code  developed  by 
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Whitlow  at  NASA  Langley  Research  Center.  Also,  several  prelimi¬ 
nary  attempts  have  been  made  to  extend  these  two-dimensional 
algorithms  to  three-dimensional  flow.  Perhaps  the  most  mature  of 
these  is  the  method  developed  by  Rizzetta  and  Borland"*-.  However, 
this  three-dimensional  code  is  very  expensive  to  run  on  a  com¬ 
puter  at  present.  All  of  these  methods  are  based  on  potential 
theory. 

In  view  of  the  computational  expense  in  obtaining  a  tran¬ 
sonic  flow  calculation,  it  is  desirable  to  increase  the  useful¬ 
ness  of  one  calculation;  to  this  end,  the  transonic  indicial 

Q  Q 

theory0' *  can  be  used.  The  indicial  theory  depends  on  the  valid¬ 
ity  of  the  principle  of  superposition.  Such  a  method  is  not 
valid  at  transonic  speeds  because  of  the  necessary  nonlinearity 
of  the  governing  equations  that  represent  the  flow  with  moving 
shock  waves.  In  certain  circumstances,  however,  unsteady  tran¬ 
sonic  flow  quantities  can  be  represented  by  a  linearized  equation 
while  the  shock  motion  is  still  taken  into  account.  Nixon9  has 
shown  that  the  transonic  indicial  method  of  Ballhaus  and 

O 

Goorjian  properly  accounts  for  shock  motion  effects,  if  only 
integrated  properties  such  as  airfoil  lift  and  moment  coeffi¬ 
cients  are  considered.  The  shock  motion  need  not  be  explicitly 
included  since  it  already  appears  in  the  calculation  of  the  indi¬ 
cial  response.  A  similar  result  applies  for  control-surface 
hinge  moments  if  the  shock  does  not  oscillate  across  the  hinge. 
When  the  shock  oscillates  across  the  hinqe,  the  strained- 
coordinate  method  of  Nixon9  can  be  used  to  treat  the  shock 
motion.  An  application  of  the  indicial  approach  is  given  in 
Reference  10. 

In  the  indicial  theory  of  Nixon,9  the  nonlinear  problem  is 
decoupled  into  two  linear  problems  in  which  superposition  is 
valid.  Since  superposition  can  be  used,  it  is  also  possible  to 
superpose  the  effects  of  various  elements  in  an  unsteady  flow. 


such  as  flap  movements  or  changes  in  angle  of  attack,  provided 
the  indicial  response  of  each  component  in  isolation  is  known. 
The  capability  to  superpose  is  one  of  the  most  powerful  features 
of  the  indicial  method. 

In  connection  with  reducing  loading  at  transonic  speeds,  it 
should  be  noted  that  a  preliminary  attempt  was  made  by  Ballhaus 
et  al .  to  study  the  effects  of  moving  leading  and  trailing  edge 
flaps  on  the  loads  on  an  airfoil. 

The  present  work  combines  the  methods  of  Kerlick  and  Nixon^0 
and  the  XTRANL  code  to  give  the  gust  response  of  an  airfoil.  In 
addition,  the  effect  of  two  control  surfaces  on  the  flow  is  in¬ 
cluded  in  the  model.  These  methods  are  augmented  by  the  use  of 
the  optimization  code  CONMIN  in  determining  the  control  surface 
deflections  to  minimize  the  additional  loading  due  to  the  gust. 
Because  the  computer  time  reguired  by  the  indicial  theory  is 
greatly  reduced,  the  optimization  code  is  economical  to  use. 

2 .  Basic  Eguations 


The  model  eguation  used  in  the  present  work  is  the  unsteady 
transonic  small  disturbance  eguation 
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where  $(x,y,t)  is  the  perturbation  velocity  potential;  x  and  y 
are  cartesian  coordinates  scaled  by  the  airfoil  chord,  c; 
and  is  the  free-stream  velocity.  The  constant  k  is  given  by 


k  =  [3  +  (y-2)  M^] 


(2) 


where  y  is  the  ratio  of  specific  heats. 
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The  appropriate  boundary  conditions  are  that  the  flow  be¬ 
haves  like  outgoing  waves  as 

2  2 

x  +  y*  ♦  -  (3) 

and  that  the  flow  is  tangent  to  the  airfoil.  This  tangency  con¬ 
dition  is  satisfied  in  a  thin  airfoil  sense  by  ensuring  that 


*y(x,±0,t)  = 


3y  (x,±0,t) 


c  3yg(x,±0,t) 


where  y  =  yg(x,±0,t)  denotes  the  ordinate  of  the  airfoil  boundary 
at  station  x  and  time  t.  The  problem  is  completed  by  ensuring 
that  the  wake,  represented  by  y  =  0,  carries  no  load.  Thus,  on 
the  wake  (and  at  the  trailing  edge  of  the  airfoil) 


Cp(x,+0,t)  =  Cp(x,-0,t) 


where  C  (x,y,t)  is  the  pressure  coefficient  given  by 


Cp(x,y,t)  -  -2(*x  +  -§-  *t) 


The  boundary  condition,  Equation  (4),  can  be  written  as 

3yc(x,±0)  3y  (x , ±0 , t ) 

♦y(x.±0.t>  =  -S-Ti—  *  -Z—x - 


e  »y3(x,±o,t) 


where  y  and  y  represent  the  steady  and  unsteady  components, 
s  s 

respectively.  For  a  control  surface  deflected  at  6 ( t ) , 

A 

yg(x,±0,t)  is  given  by 


yg(x,±0,t)  =  -B(t)  (x-xh)  H(x  -  xh) 


-5- 


where  is  the  location  of  the  control  surface  hinge  and  H(  )  is 
the  step  function. 


In  the  present  analysis,  a  semi-infinite  gust  is  modeled  by 
the  boundary  condition 


<fry(x,±0,t)  = 


3y  ( x ,±0 ) 


+  vG(x,t) 


Thus, 


vG(x,t)  =  |  vG|  ,  x  <  UJt  -  tQ) 

=  0  ,  x  >  UJt  -  tQ) 


where  |vG|  is  the  magnitude  of  the  gust  velocity;  the  gust  starts 

at  t  =  t0  and  travels  with  the  free-stream  velocity  U^.  The 
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gust  model  is  similar  to  that  of  McCroskey  and  Goorjian* 

3.  Piston  Theory  Limit  for  a  Step  Change  in  Normal  Velocit 


In  the  early  phase  of  subsonic  and  supersonic  flutter  analy¬ 
sis,  great  use  was  made  of  the  piston  theory  limit.  Briefly,  the 
supposition  was  that  if  the  airfoil  oscillated  at  very  high  fre¬ 
quencies,  the  velocity  of  the  airfoil  surface  was  an  order  of 
magnitude  greater  than  the  streamwise  velocity.  As  a  conse¬ 
quence,  the  linearized  version  of  Equation  (1)  becomes 


<t>  =  M2 

yy  -  02  ytt 


The  solution  is 


For  y  >  0,  the  positive  sign  is  taken;  for  y  <  0,  the  nega¬ 
tive  sign  is  used.  This  is  to  ensure  outgoing  waves.  The  pres¬ 
sure  coefficient  in  the  approximation  of  Equation  (11)  is  given 
by 


Cp(x,±y,t)  =  -2 


U 


4>t(x,±y,t ) 


(13) 


From  Equation  (12), 


♦y(x,y,t)  -  *■( y  +  -jf-  t) 


(14) 


For  a  step  change  in  the  normal  velocity  component  of  -a,  Equa¬ 
tion  (14)  gives 


f  '  (0)  =  -« 


(15) 


Combining  Eguations  (13)  and  (15)  gives  the  result 

C  (x ,±0,0)  -  ?  (16) 

This  is  the  initial  pressure  coefficient  experienced  by  the  air¬ 
foil  due  to  a  step  change  in  angle  of  attack. 

In  the  transonic  equation.  Equation  (1),  there  is  no  mecha¬ 
nism  to  invalidate  the  basic  premise  of  piston  theory;  namely, 
that  streamwise  velocities  are  small  relative  to  the  normal  velo¬ 
city.  Hence,  Equation  (16)  is  valid  for  transonic  flow.  Fur¬ 
thermore,  since  streamwise  velocities  contain  the  shock  mechanism 
for  Equation  (1),  the  initial  movement  of  shock  does  not  enter 
Equation  (11).  Therefore,  it  is  assumed  that  the  shock  movement 
does  not  start  with  the  discontinuity  that  characterizes  the 
pressure  coefficient. 


r»  '--v  i  i  m.  r_  i  .  .  .  .  »  .  .  .  .  .  ■  -  -  * 


4.  The  Transonic  Indicial  Method 


In  linear  models  of  subsonic  or  supersonic  flow,  an  alterna¬ 
tive  to  expanding  the  governing  eguation  in  a  Fourier  series  is 
to  use  an  indicial  function.  An  indicial  function  analysis  de¬ 
pends  on  the  governing  eguation  being  linear,  which  is  not  the 
case  for  transonic  flow.  However,  if  the  amplitude  of  motion  is 
small  and  if  shock  waves  are  not  generated  or  destroyed  during 
the  motion,  a  linear  eguation  similar  to  that  used  in  subsonic 
flow  can  be  derived  to  model  the  unsteady  motion.  This  is  the 
basis  for  the  work  of  Kerlick  and  Nixon. The  following  briefly 
outlines  their  work: 


Consider  the  change  in  the  pressure  coefficient  Cp(x,y,t) 
due  to  an  infinitesimal  change  in  some  parameter  e(x)  at  some 
time  t.  If  the  pressure  distribution  varies  continuously 
with  e,  a  Taylor's  series  expansion  gives 


ACp(x,y,t) 


Cp  (x ,y , t , x ) 

e 


de  (  t  ) 
dt 


At 


higher-order 

terms 


(17) 


where  t  is  the  time  and  C  (x,y,t,x)  is  the  rate  of  change  of 
C  (x,y,t)  with  e  at  some  time  t;  AC  (x,y,t)  is  the  change  in 

M  p 

pressure  distribution.  Neglecting  the  higher  order  terms,  the 
total  effect  of  all  such  steps  up  until  time  t  is  then 

AC  (x,y,t)  =  e(0)C  (x,y,t,0)  +  f C  (x,y,t,x)  d*(t)dx  (18) 

P  Pe  0  Pe  ax 

If  it  can  be  assumed  that  the  behavior  of  Cp  with  e  is  linear, 

C  (x,y,t,x)  can  be  represented  by  its  value  at  x  =  0, 

*  e 

provided  that  the  time  t  is  taken  relative  to  t.  Thus, 


C„  (x,y,t,T)  =  C  (x,y,t-x,0) 


(19) 


Equation  (18)  then  becomes 


AC  (x,y,t)  =  e(0)C  (x,y,t,0)  +  /  C  (x,y,t-T,0) 


de  ( t ) . 


(20) 


By  a  simple  change  of  variable,  Equation  (20)  can  then  be 
written  as 

AC  (x,y,t)  =  e ( 0 )C  (x,y,t)  +  /  C  (x,y,T)  — ■  -T—  dt  (21) 

P  Pe  0  pe  QT 

where  the  functional  form  of  C  (x,y,t,0)  has  been  contracted  to 

Pe 

C  (x,y,t)  for  convenience  in  presentation.  Equation  (21)  is  the 
Pe 

form  of  a  typical  relation  in  the  indicial  method.  It  is  assumed 

that  C  (x,y,t)  can  be  obtained  either  analytically  or  numeri 
Pe 

cally.  In  the  foregoing  analysis,  it  has  been  assumed  that  the 
variation  of  Cp(x,y,t)  with  e  is  linear.  For  continuous  flows 
governed  by  a  nonlinear  equation,  this  assumption  may  well  be 
adequate  for  small  values  of  e.  However,  as  pointed  out  by 
Nixon,9  if  there  are  shock  waves  that  move  durinq  the  motion,  the 
assumption  of  linear  variation  is  totally  inadequate  in  the 
region  bounded  by  the  extremities  of  the  shock  motion.  This  is 
because  the  pressures  in  this  region  can  jump  from  a  pre-shock 
value  to  a  post-shock  value,  or  vice  versa,  as  the  shock  tra¬ 
verses  this  region.  This  difficulty  can  be  overcome  by  the  use 
of  a  strained  coordinate  system  in  which  the  shock  remains  at  the 
same  location. 

Let  the  strained  coordinate  system  be  given  by  (x',y,t), 
which  is  related  to  the  physical  coordinate  system  (x,y,t)  by 


x  =  x'  +  fix  (t)x, (x' ) 
s  l 


where  fix  (t)  is  the  change  in  shock  location  at  some  time  t,  and 
s 

xj(x')  is  a  known  straining  function.  In  Equation  (22),  it  is 
assumed  that  the  shock  waves  are  normal  to  the  free  stream. 


Following  the  ideas  of  Nixon,  Cp(x,y,t)  is  given  by 


C  (x,y,t)  =  C  (x',y)  r 1  -  6x  ( t )x,  (x’)l  +  C  (x',y,t) 

P  P0  l  8  XX'  J  P1 


(23) 


where  Cp^(x',y)  and  Cp  (x',y,t)  are  the  values  for  the  mean 


steady  state  pressure  and  the  unsteady  perturbation, 
respectively.  The  coordinates  x'  and  x  are  related  by  Equation 
(22). 


In  Equation  (23),  6xe(t)  and  C  (x’,y,t)  are  linearly  depen 

S  P-t 


dent  on  the  parameter  e(t)  and,  therefore,  can  be  treated  by  the 
indicial  method.  The  nonlinear  effect  appears  implicitly  through 
the  transformation  from  the  strained  coordinates  to  the  physical 
coordinates.  Thus, 


6x  ( t )  =  6xe  (t)  e  ( 0)  +  /  5x  (t)  lSit_ILJLldT 
s  Se  0  se 


dt 


(24) 


c  (x',y,t)  =  c  (x',y,t)  e ( 0)  +  /  c  (x ' ,y , t )  T)dt  (25) 

P1  Pe  0  Pe 


dx 


where  6x  (t)  and  C_  (x',y,t)  are  the  indicial  responses  of  6x  (t) 
sp  s 

E  E 


and  Cp^(x',y,t),  respectively.  The  indicial  responses  can  be 


obtained  by  computing  the  transient  behavior  of  the  shock  motion 
and  the  pressure  coefficient  due  to  a  step  change  in  e.  This 
transonic  calculation  must  only  be  performed  once  for  a  range  of 
functions  e(t). 


Equations  (24)  and  (25)  are  linear;  thus,  the  principle  of 
superposition  can  be  used.  In  the  present  work,  the  unsteady 
flows  due  to  a  gust  and  two  control  devices  are  superposed.  In 
this  case. 


T 


(26) 


3  (  t  de.(t-r) 

«xs(t)  mt  Uxs  (tJejCO)  +  /  6xg  (t)  — ^ - d 


i  =  1 


0  e; 


C  (x',y,t)  =  l 


i  =  1 


C  (x' ,y,t)ei(0) 
ei 


t  de .  (t--r) 

+  /  C  (x',y,T)  — Vr - dr 


0  Pe; 


dt 


(27) 


where  the  subscript  i  denotes  a  flow  perturbation  due  to  the  ith 
device  (either  controls  or  the  gust). 


5.  Curve  Fits  of  Indicial  Data 


It  is  considerably  helpful  in  the  numerical  evaluations  of 
Equations  (26)  and  (27)  if  the  indicial  responses  can  be  repre¬ 
sented  by  analytic  functions  since  the  convolution  integrals  can 
then  be  evaluated  analytically.  In  the  present  work,  the  indi¬ 
cial  responses  for  both  the  pressure  coefficient  and  the  shock 
movement  are  approximated  as 

f(x,t)  =  a  +  b  exp(-at)  +  c  exp(-Bt)  (28) 

where  f(x,t)  is  a  general  indicial  response  and  a,  b,  c,  a, 
and  B  are  functions  of  x,  in  general.  The  reason  for  including  a 
second  exponential  term  in  Equation  (28)  is  to  allow  for  an  in¬ 
flection  point  on  the  indicial  response  curve.  At  t  *  0,  f(x,0) 
is  represented  by  its  piston  theory  limit  fp(x).  Thus, 

f p ( x )  =  a  +  b  +  c  (29) 

As  t  ♦  f(x,t)  must  approach  its  steady  state  value  f(x,«) 
corresponding  to  the  step  change  in  e.  Thus, 


I 


f(x,-)  =  a  (30) 

Combining  Equations  (28),  (29),  and  (30)  gives 

f(x,t)  *  f(x,«)  II  -  exp(-ot)]  +  £p(x)  exp (-at ) 

+  c  lexp(-Bt)  -  exp(-at)]  (31) 

Apart  from  the  initial  jump  predicted  by  piston  theory,  the 
indicial  response  is  a  smooth  curve.  Consequently,  [£(x,t)  - 
fp(x)]/f (x,«)  can  be  approximated  by 

[f ( x , t )  -  fp(x)]/f (x,»)  «  [l  -  fp(x)/f (x,«)]  [l  -  exp(-at)] 


+  fY x , » )  [exp(-Bt)  -  exp(-at)]  (32) 

The  actual  values  of  £(x,t)  are  generated  by  numerical  data. 
The  coefficients  c,  a,  and  6  are  found  by  using  Equation  (32)  in 
conjunction  with  the  computer  program  CONMIN,  which  minimizes  the 
objective  function 

[ f  ( x , t )  -  fp(x,0)] /f  (x,«)  -  [l  -  f  (x)/f  (x,«)]  [l  -  exp(-at)] 

-  c/f(x,»)  [exp(-Bt)  -  exp (-at )] 
over  the  computed  range  of  t.  For  the  pressure  coefficient 


fp(x) 


—  +2a 
M 


and  for  the  shock  motion 


fP- 0 
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The  program  for  the  curve  fit  was  based  on  an  early  version  pro¬ 
vided  by  Dr.  Samuel  C.  McIntosh  of  McIntosh  Structural  Dynamics, 
Inc . 

6 .  Model  of  Gust  Alleviation 

A  model  problem  was  developed  to  test  the  applicability  of 
the  indicial  method  to  model  gust  alleviation.  The  problem  con¬ 
sists  of  a  "top  hat"  gust  that  crosses  the  leading  edge  at  time 
tjQ  and  two  control  surfaces.  For  generality  of  computation,  the 
sharp  edge  of  the  gust  is  replaced  by  a  continuous  variation 
similar  to  that  of  the  control  surfaces  as  given  below.  The 
characteristic  times  of  the  control  surface  are  t10,  tj^,  tj2, 
and  t13. 

At  t  =  t20,  the  control  surface  moves  according  to  the  law 
a  =  aQ  [1  -  cos  ~  t2Q)I/2,  t2Q  <  t  <  t21  (35) 

where  a>21  is  given  by 


“21  =  *^*21  “  fc20* 

At  time  t2i  the  control  stops.  At  time  t22,  the  control  moves 
according  to  the  law 


a  ■  aQ  (cos  «>22(t  -  t22)  +  l]/2,  t22  <  t  <  t23  (37) 


where 


The  second  control  surface  moves  in  an  analogous  manner  on  the 
first  control. 

The  control  motion  described  by  Equation  (35)  has  zero  angu¬ 
lar  velocity  at  t  =  t2o  and  t2^,  which  seems  to  be  a  realistic 
assumption.  A  similar  behavior  occurs  at  t22  and  t2j . 

In  the  indicial  formulation,  a  typical  quantity  is  given  by 

f(t)  =  f  ft)  e  ( 0)  +  /%  (t)  de^~T-)dT  (39) 

e  q  e  at 

where  f£(t)  is  the  indicial  response.  For  the  control  laws  de¬ 
scribed,  f(t)  is  given  by 

A 

°o  fcl 

f(t)  =  — 2  u>  /  f£(t)  sin  u(t  -  t  +  tp)  dt  (40) 

^0 

where  the  generic  control  law 

[1  -  cos  w(t  -  tQ)]/2  (41) 


is  used. 

If  the  control  surface  movement  starts  at  tQ  and  finishes  at 

tl' 

~CdT~~'  =  0  for  t  -  t  +  tQ  <  tQ  (42) 


or,  for  a  nonzero  value  of  the  integrand, 

t+tQ-tj<T<t  (43) 


Also,  for  the  general  limits  on  the  integral. 


I 


J. 


Hence,  the  limits  are 


tQ  =  max(tQ,  t  +  tQ  -  tx  ] 

(45) 

* 

fcl  =  fc 

7 .  Control  Optimization 

Calculations  using  the  indicial  response  require  very  little 
computer  time.  Thus,  it  is  economical  to  couple  the  gust 
response  and  the  control  model  to  an  optimizer  to  determine  the 
best  control  deflections  and  to  minimize  the  objective  function 
F,  where 

F  =  AC^  +  u  AC*  (46) 


There,  ACL  and  ACM  are  the  incremental  lift  and  moment  coeffi¬ 


cients,  and  w  is  a  specified  weighting  function.  The  moments  are 
about  the  quarter  chord.  The  variables  used  in  the  optimization 


are  the  times  t^k  (j  =  2,3;  k  =  0,1, 2, 3)  described  in  the  pre¬ 


vious  section  and  the  amplitudes  of  the  control  surface  deflec¬ 
tions.  Constraints  on  tjk  are  enforced;  namely,  that  the  control 
surface  cannot  operate  until  the  gust  has  passed  by  0.1  time 
units  and  that  the  finish  time  of  a  control  deflection  cannot 
occur  before  the  start  time  of  the  deflection. 


The  objective  of  this  exercise  is  to  devise  simple  control 
systems  to  counter  the  gust.  The  computer  code  CONMIN  is  used 
for  the  optimization. 


8.  Gust  Response 


In  linear  subsonic  flow,  the  gust  response  is  a  function 
only  of  the  free-stream  Mach  number  and  the  gust  character¬ 
istics.  In  nonlinear  transonic  flow,  the  response  is  also  a 
function  of  the  steady  flow  over  the  airfoil  and,  hence,  the 
airfoil  geometry.  This  coupling  manifests  itself  in  two  ways. 
First,  the  asymptotic  lift  of  the  wing  due  to  a  step  in  angle  of 
attack,  a,  can  be  much  greater  than  the  classic  value 
of  2jra/8,  where 


Second,  the  characteristic  time  of  the  response,  that  is  the 
transient  behavior,  is  different  from  that  in  subsonic  flow.  For 
an  airfoil  to  have  a  low  gust  response,  the  asymptotic  lift 
increment  of  CL  due  to  the  gust  should  be  small.  Also,  the  char¬ 
acteristic  time  should  be  large  to  allow  for  a  control  to  oper¬ 
ate. 


For  a  gust  alleviation  device,  the  incremental  lift,  ACL, 
should  be  large  so  that  the  control  is  effective;  the  charac¬ 
teristic  time  should  be  small  so  that  the  effect  of  the  control 
occurs  as  early  as  possible. 


In  transonic  flow  there  is  a  nonlinear  coupling  between  the 
steady  and  unsteady  flows  that  possibly  could  be  used  to  some 
advantage.  Consider  that  the  unsteady  lift  coefficient  CL  is 
represented  as  follows: 


For  a  gust 


C 


L 


f 

00 


r 


-c(t  +  tft)l 


/  a  n  \ 


where  f  is  the  value  of  CL  as  t+<»,  c  is  a  constant,  and  the  gust 
starts  at  t  =  -tg. 

For  a  control  surface 

C.  =  f  +  (f  -  f  )  e~ct  (48) 

L,  oo  p  oo 

where  f^  is  the  value  of  CL  at  t-*-»,  and  fp  is  the  value  for  pis¬ 
ton  theory  at  t  =  0. 

Let  the  control  law  be 


o  =  — j  (1  ~  cos  ojt ) 


0  <  t  <  t. 


t  >  t. 


Olt^  =  IT 

Using  Equations  (47),  (48),  and  (49)  in  the  indicial  equations  of 
section  4  gives  the  total  unsteady  lift  as 

-ctf(l  +  eY")  (F  -  1 


where 


F 


g  e  0  (1+Y  ) 


=  1  +  2  - 

(1  +  eYir) 


or,  using  Eguation  (52), 


F 


e-ct0  (1+Y2) 
(1  +  e Y?r ) 


(53) 


For  linear  subsonic  flow  at  a  fixed  Mach  number,  f^,  f  ,  and  c 
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are  independent  of  airfoil  section?  u>  is  specified.  Hence,  there 
is  no  mechanism  by  which  Eguation  (53)  can  be  satisfied  in  gen¬ 
eral.  However,  for  transonic  flow  the  guantities  f  and  c  are 

CO 

dependent  on  airfoil  geometry,  and  it  is  possible  that  a  certain 
design  of  airfoil  and  control  surface  could  satisfy  Equation  (52) 
leading  to  an  effective  gust  alleviation  device. 


9.  Results 


The  computer  code  XTRANL  is  used  to  compute  the  necessary 
indicial  responses  and  to  compute  fully  nonlinear  cases  for  com¬ 
parison  . 

In  the  following  cases,  the  gust  is  of  a  "top  hat"  geometry, 
starting  at  t  =  0  and  ending  at  t  =  50.  There  are  two  control 
surfaces,  each  of  10-percent  chord  extent.  One  control  is  at  the 
leading  edge,  the  other  is  at  the  trailing  edge.  In  the  test 
cases,  the  controls  operate  at  the  time  shown  in  the  following 
table.  For  the  optimization  cases,  the  gust  amplitude  is  1.5 
degrees . 


Table  1.  Values  of  Trj 


IS 

0 

1 

2 

3 

0.0 

4 . 36xl0-5 

50.0 

50+4 . 36xl0~5 

2 

0.0 

22.0 

60.0 

82.0 

3 

0.0 

22.0 

60.0 

82.0 

Note:  r  =  1,  Gust 

r  =  2,  Trailing  edge  control 
r  =  3,  Leading  edge  control 


The  steady  pressure  distribution  around  a  NACA  0012  airfoil 
at  a  =2.89  degrees  and  M  =  0.601  is  shown  in  Figure  1.  This 

iti  ® 

is  a  subcritical  example.  In  Figure  2,  the  prediction  of  the 
lift  for  this  steady  configuration  is  shown  with  a  "top  hat"  gust 
with  amplitude  eguivalent  to  0.25  degrees  and  two  control  sur¬ 
faces  deflected  at  0.25  degrees.  The  agreement  with  the  direct 
result  is  guite  good. 

In  Figures  3  and  4,  the  incremental  lift  and  moment  coeffi¬ 
cients  for  the  gust  are  shown.  The  function  F  in  Eguation  (46) 
is  minimized  through  the  overall  time,  and  the  optimized  incre¬ 
mental  lift  and  moment  coefficients  are  also  shown.  In  the  opti¬ 
mization,  the  magnitude  of  the  lift  coefficient  was  much  greater 
than  that  of  the  pitching  moment,  and  the  function  F  represents 
essentially  only  the  lift  increment.  Accordingly,  to  make  both 
lift  and  moment  terms  of  comparable  size,  the  incremental  lift 
and  moment  were  scaled  with  respect  to  the  peak  lift  and  moment 
increments  for  the  gust  alone.  These  results  are  also  shown  in 
Figures  3  and  4.  It  can  be  seen  that  in  the  present  simple 
analysis  it  is  difficult  for  the  two  controls  to  alleviate  both 
the  lift  and  moment,  although  it  is  possible  to  control  the  lift 
itself  guite  well. 


A  similar  series  of  data  is  shown  in  Figures  5  through  8  for 
a  NACA  64A006  airfoil  at  a  =  1.0  degree  and  M#  =  0.825.  Again, 
the  agreement  between  the  direct  result  and  the  present  result  is 
good.  In  the  optimization  cases,  similar  conditions  to  those 
noted  for  the  previous  example  apply. 

A  third  transonic  case  is  shown  in  Figures  9  through  12. 

The  airfoil  is  an  MBB-A3  section  at  a  =  1.5  degrees  and  M  = 
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0.7.  The  agreement  between  the  direct  result  and  the  present 
result  is  good;  the  same  general  conclusions  regarding  optimiza¬ 
tion  noted  earlier  apply. 

The  fact  that  the  gust  loading  of  the  first  transonic 
example  can  be  alleviated  more  than  the  subsonic  example  would 
reinforce  the  suggestion  made  in  section  8. 

10.  Possible  Extensions  of  Work 


The  main  advantage  of  the  present  technique  is  that  the 
computational  effort  is  minimal,  thus  allowing  the  use  of  an 
optimization  routine  to  determine  the  best  control  movements. 
Provided  the  number  of  shock  waves  in  the  flow  does  not  change, 
the  method  is  very  useful.  If  the  number  of  shock  waves  does 
change  during  the  unsteady  motion,  the  present  theory  is  invalid; 
however,  it  may  be  possible  to  extend  the  theory  to  treat  this 
aspect.  Such  an  extension  for  steady  flow  is  given  in  Reference 
13.  There  is  no  difference  in  principle  in  the  indicial  theory 
between  two  and  three  dimensions,  and  a  logical  extension  of  the 
present  work  would  be  to  three-dimensional  flows. 


Several  aspects  of  the  existing  two-dimensional  procedure 
need  improving.  One  aspect  concerns  the  problems  that  can  arise 
when  the  supersonic  bubble  is  small  and  near  the  leading  edge. 
The  present  method  uses  only  one  point  straining  (at  the  shock 


wave).  It  has  been  found  in  earlier  work  on  the  steady  flow 
problem  that  the  points  characterizing  large  gradients  also  need 
to  be  strained.  Such  a  modification  is  a  straightforward  task. 
A  second  aspect  is  that  the  optimization  code  CONMIN  may  not  be 
the  best  code  to  use;  further  study  of  the  possibilities  of  using 
a  different  optimization  code  is  desirable. 


Because  of  the  computational  speed  of  the  present  method,  it 
can  be  used  to  investigate  different  methods  of  gust  alleviation, 
including  radical  concepts.  The  main  restriction  on  the  present 
work  is  that  the  aerodynamics  of  such  devices  must  be  capable  of 
being  modeled  by  potential  theory.  This  restriction  is  because, 
at  present,  unsteady  transonic  flows  are  not  routinely  calculated 
by  the  Euler  or  Navier-Stokes  equations.  If  such  simulations 
were  available,  the  indicial  method  would  be  applicable.  It 
should  be  noted  that  a  Navier-Stokes  simulation  is  similar  to 
experimental  data  in  principle. 

11 .  Conclusions 

A  method  for  predicting  the  gust  loading  on  an  airfoil  mov¬ 
ing  at  transonic  speeds  based  on  the  indicial  theory  has  been 
developed  that  is  computationally  efficient.  The  method  allows 
easy  incorporation  of  gust  alleviation  devices  and,  when  linked 
to  an  optimization  method,  gives  the  best  mode  of  operation  of 
these  control  devices.  The  method  should  be  applicable  to  three- 
dimensional  flows  where  a  considerable  savings  in  computer  re¬ 
sources  over  a  direct  simulation  of  the  unsteady  flow  could  be 
obtained.  The  method  is  semi-analytic  in  nature,  which  allows 
greater  insight  into  methods  of  gust  alleviation  than  can  be 
obtained  from  a  purely  numerical  study.  Further  study  of  the 
equations  should  lead  to  new  ideas  for  alleviating  the  effect  of 
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